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The hierarchy of Green functions for (quasi) degenerate systems, presented in cond-mat/0308058, 
is calculated in detail for the case of a system with closed shells plus a single electron in a two-fold 
degenerate level. The complete hierarchy is derived explicitly, in terms of Green functions and of 
Feynman diagrams. 

I. INTRODUCTION 

In reference 1 , hereafter referred to as I, a many-body theory of degenerate or quasidegenerate systems was presented, 
using quantum group methods. In particular, a hierarchy of Green function was derived and formulated in terms of 
reduced coproducts. In this paper, the general formula is written explicitly in terms of Green functions and Feynman 
diagrams. The case we consider is that of a system with closed shells and a two-fold degenerate level occupied by a 
single electron. This does not look like a very interesting example, but we shall see that the hierarchy involves n-body 
Green functions with n up to 5. In the general case of a M-fold degenerate level, the Green functions involved go up 
to n = 3M — 1. In practice, such complicated Green functions are not usable, and the hierarchy is stopped after the 
first terms. The hierarchy we give here contains the first terms of the general hierarchy. 

Quantum group concepts are not (yet) familiar in solid-state physics, so the pace will be slow and the calculations 
will be presented in full detail. 

II. CALCULATION OF W x p 

The simplest case is one electron (N = 1) in a twofold degenerate level (M = 2). The two degenerate levels are 
denoted by 1 and 2. The density matrix is 

( Pn P12 

\ P21 P22 

with tr(p) pn + p 22 = 1. 

The moment generating function p(a,a) is defined in equation (28) of I. In our case, N = 1 and p(a, a) = 
po(a, a) + pi(a, a). According to equation (29) of I (with N = 1, i\ = 1 or 2, ji = 1 or 2) 

Pl(a,Ot) = p\\OL\OL\ + p\2&\OL2 + P2\&20L\ + P220L20L2- 

Using now equation (30) of I with k = 1 we obtain 

po(a,a) tr(p) = pn +p 2 2 = 1- 
According to equation (32) of I, the cumulant generating function p c (a,a) is obtained from 

log(p(o;,a)) = log (l + pi(a,a)) = pi(a,a) + p c (a,a), 

where 

p (ci Ot)^ 1 

p c (a,a) = p c 2 (a,a) = ^ = - det (p)a 1 a 1 a 2 a 2 = -{pnp22 - pi2P2i)«i«i«2«2. 

According to equation (34) of I, this enables us to calculate W®. But we saw in section 7.3 of I that Wl = W®, thus 



with 

C 2 2 



G° p (x,y) = G° (x,y) + i(j2 

(x)u m (y) + ^2^2pi k u k (x)ui(y)^)M, 



k=i i=i 
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where 

M 



1 -1 
-1 1 



The expression for the 2x2 matrix Gq(x, y) is given in equation (15) of I. 
For later convenience, we give the identities 

S^ 1 o^i o^i cy, 2 (^2 

cc'(ui(x)a 2 - u 2 (x)ai)(ui(y)a 2 - U' 2 (y)ai). 



5rj e >(y)5r] e (x) 

5 A p c 2 (a,a) 



-e 1 €2e' 1 e' 2 r(x 1 ,X2,yi,y2), (2) 



Srj €1 (x!)Srf e2 (x 2 )Sr) € f i (yi)5r) € > 2 (2/2) 
where 

r(x 1 ,X2,yi,y2) = - det(p)(u 1 (xi)u2(x 2 ) - u 2 (x 1 )u 1 (x 2 ))(ui(yi)u2(y2) ~ ^2(2/1)^1(2/2))- 

III. HIERARCHY 

For one electron in a two-fold degenerate state, equation (43) of I stops at n = 3 because M = 2 and we have 



To obtain the Green function, we put /? = fj e (x), we make a functional derivative with respect to r/ e /(?/) and we set 
all fermion sources to zero (in particular, this implies Z p = 1). We use also the fact that a functional derivative of 
Z p or Wp with respect to an odd number of fermion sources is zero if the sources are set to zero. This enables us to 

calculate the signs ( — 1) . 

52 Z " 52W ? ^(D ^ VP z) i^(D SW })(D 5Z " ) 



A. Coproduct 



The non-relativistic interacting Hamiltonian is 
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| ff int (t)dt = y J Vi(i,r)^(i,r')^^V(*,r')V'(i,r)dMrdr'. (4) 

To simplify the notation, we define #1 = (t, r), = (£', r') and v(xi — x^) = 5(£ — t')e 2 /|r — r'|, so that 

J H int (t)dt = X - J ^(^-^^^^^(^^(^^(xOd^dx;. 

To define the interacting operator D (i.e. the functional differential form of J H int (t)dt), we make the substitutions 
ip(xi) — > —iS / Sfj(xi) and ^(xi) — > iS/Srj(xi) and we use the simplified notation 

tfi = g = g ^ = g £ ^ 



^±(xi)' ± Srj±(x 1 )' ± ^±K)' ± ^?±0i)' 



This gives us D = D+ — D_ with 

D± = ^ J v(x 1 - x^S^Si^S^S^dxidx^. (5) 

Thus 

D = \ J d^idxi^Oi - - Si_s v Jvji^, 

To further simplify the notation, we replace 5\ ± by Si and x\ by 1 and by V . The interacting operator becomes 

D = - [ dldl'v(l-l')^2±5 1 S v Si,5 1 . 

To calculate the coproduct of D, we start from the coproduct of the basic functional derivatives 

ASi = Si 1 + 1 5i, AJi/ = (Ji/ 1 + 1 A5i = 5i 1 + 1 (8) <*i, A^, = 1 + 1 

Then we use equation (21) of I to calculate the coproduct of SiSy and SySy 

ASiSy = SiSy 1 + Si Sy - Sy S\ + 1 ASySi = 1 + (Sp <Si - <Si tfp + 1 c^i. 

We use equation (21) of I again to obtain 

A<Mi'£i'£i = ^l^i^i/^i 1 + 10 + SiSySit 5i - 5i/ + + 5i£i'<*i 5i/ 

— SiSy <5j + SiSi + 5i SySySi — SySySi Si + — <5i/<5i tfitfp 

<Mi/<5i + + 5i/ - Si <Mi/£i/. (6) 

The reduced coproduct with respect to D is obtained by its definition A'D = AD - Therefore, 
A'D = dldl^(l - l , )^2±(s i SySy Si - SiSySi Sy + frfr; SySi + ^l^i/^i Sy 

—SiSy + SiSi (5i/ (5^/ + 5i SySySi — SySySi Si + ^/^p — 

+ SiSy + 5i/ - 5i (7) 

Notice that, in A'D = X^^a') ®-^(2')> eacn term -D ( i') or -D( 2 ') contains between one and three functional derivatives. 

B. Expansion of equation (3) 

Now we are going to examine each term of equation (3). We define the n-body Green function as 

S 2n Z 

Ge 1 ...e„e' 1 ...e' n (xi, . . . , X n , y[, . . . ,y' n ) = (-<)» — ^ ■ ■ ■ — j- ^ ■ ■ ■ — ^ . (8) 

The right hand side is taken at zero external sources. To simplify the notation, we write i for X{ and z' for we 
consider that goes with the variable a^, e- with ^, so the e are now implicit and we have 

G(l,...,n,l',...,n') = (-i) n S-i...SnS v ...Sn>Z p . 

With this notation, the left hand side of equation (3) is now 

s 2 z P s 2 z p 

- -iG(x,y). 



Srj € >(y)Sr) € (x) Sr] e (x)Sri e >(y) 

More explicitly G(x,y) = G €€ /(x,y). In the following we shall calculate in detail all the terms of the right hand side 
of equation (3). 
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C. First term of equation (3) 

The first term of (3) is very simple. From equation (1), we know that Wp is the sum of a term of degree two and 
a term of degree four in the external sources. The first term in the right hand side of equation (3) contains only two 
functional derivatives. Thus, at zero sources, the second term of Wp does not contribute. By definition of functional 
derivatives, the first term gives us 



Pw 1 

Sr) e f(y)Sr) e (x) 9 



D. Second term of equation (3) 

The second term of equation (3) is 

-iY,{D (ll) 5 y 5 s Wl){D (2l) Z p ), 

where we wrote S y for 5/5r] e >(y) and 5$ for S/Sfj e (x). According to equations (7), each term D {1 > ) contains at least one 
functional derivative. Therefore, the first term of Wp 1 does not contribute because it is of degree 2. The second term 
of Wp is of degree four, and it contains two sources r] and two sources fj. The functional derivatives are evaluated at 
zero sources, thus D (1 ^ must be the product of a derivative with respect to a source r] and a derivative with respect 
to a source rj. According to equation (7), there are four terms with this property in A'D: 

\ /dldiyi - i')I>(- ^ ® w + ® *i>h. + ® S lSi - ^ ® Mi) 

The second term of equation (3) becomes 

- l - J dldl'i;(l - l')X)±(- (Si8i>6yS s W})(8v6iZ p ) + (S^SySsW^iSySyZp) + {5 v 5 v 5 y 5 s Wl) {S^Zp) 

-(SvSidySxW^^SyZp)). (9) 

Let us now consider the first term of (9). According to (8), SyS^Zp = — iG(l, 1'). According to (2), 

S^SyS^Wp 1 = S^- v S 1 S y Wp 1 = -ee f r(x,l , ,l,y). 

For the second term of (9) we have 5y5yZ p = — zG(l', 1') and 

hk&y&xWl = 5x5-^^1 = -ee'r{x,l,l,y). 

Interchanging the variables 1 and 1/, we see that the third term of (9) is equal to the first one, and the fourth term 
is equal to the second one. Finally, the second term of (3) becomes the sum of 



and 



-ee' f dldlVl-l')r(^,l',l,y)^±G(l,l'), 
■' ± 

' J dldl't;(l - l')r(x, 1, 1, y) ±G(1', I')- 



E. Third term of equation (3) 



The third term is 



-iY,{P«f*Wl){D^8 y Z p ). 
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For this term, two kinds of D {1 > ) ® D {2 > ) intervene: with one or three products in D (1 > y For the first kind of terms, 
the relevant coproducts are 6\ 6y6y6\ — 6y 6i6y6\ and the result is 

- % - J dldiyi - I') ±((Si6xW}) (6y6y6 l 6 y Zp) - (SvSzWfiiSiSiiSiSyZ,,)) 

= -i [ dldl , v(l-l , )^2±(5 1 5^)(S v 5 v 5 1 S y Z p ). 

We have interchanged the variables 1 and 1' to obtain the last line. We use now 6\6xWp = —iG° p (x, 1) and 

6y6y6\6 y Z p = — 5\5y5y5yZ p = G(l, 1', 1', y), 

to get 

- f dldiy 1 - 1') ±G° p (x, 1)G(1, 1', 1', y). 

For the second kind of terms, the relevant coproducts are 6\6y6y ® Si — 6i6y6i <g> 6y. Again, interchanging 1 and V 
shows that these two terms give the same contribution and the result is 

i [ dldl , v(l-l , )^2±(5 1 5 v 5 1 S^)(5 v S y Z p ). 

We calculate 

6 1 6y6 1 6 x W p 1 = -5*5 l 5 1 5 1 ,W} = ±er(x,l,l,l'), 
and 8y6 y Z p = iG(¥, y). The contribution of the second kind of terms becomes 



• j dldl't;(l - \')r{x, 1, 1, 1') V)- 



F. Fourth term of equation (3) 

The fourth term of equation (3) is 

-\J2( d I'^ w p)( d U z p )- 

Here we meet something new because we must use the reduced coproduct A'D 2 . According to the general formula 
(39)of I, this reduced coproduct is obtained by multiplying A'D by itself. The spacetime variables of the first A'D will 
still be denoted by 1 and 1', while the spacetime variables of the second A'D will be denoted by 2 and 2'. Moreover, 
the ± variable of the second A'D will be written ±' . Each term D {1 > ) of A'D is at least of degree 1, so each term 
of D 2 ±/) is at least of degree 2. Therefore, the only term of W p that contributes to D 2 1 , ) 6 y 6xWp is p c (a,a). We use 
again the fact that p c (a, a) is a product of two sources r] and two sources fj to deduce that D 2 lf) must be the product 
of a functional derivative with respect to r] and a functional derivative with respect to fj. The eight terms D 2 lf) that 
satisfy this condition are 5i8 2l SiSy, 61*62, 6y8 2 ' 1 ^1^2? ^1^2? 6162', 61,82 and SySy. This gives us eight terms which 
are all identical if we interchange 1 with 1' or 2 with 2' or (1,1') with (2,2'). Thus, the fourth term of equation (3) 
becomes 

- / dldl / d2d2 / v(l - 1>(2 - 2') ± ±' (SiS2SyS^Wp 1 )(SySyS- 1 S 2 S 2f S^Zp). 
We have just to calculate 

SyS^SyS^Wp 1 = S^SySyWp 1 = -ee' ± ±'r(x, 2, 1, y), 

and 

SiSy 616262' SyZp = 621 6\6y 6262' 6yZ p = -iG(2 / , 1, 1', 2, 2', l'). 
The final expression for the fourth term is then 

-166 1 j dldl / d2d2 / v(l - 1>(2 - 2')r(x, 2, 1, y) ^ G(2 / , 1, V , 2, 2', l'). 
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G. Fifth term of equation (3) 

The fifth term of equation (3) is 

-\T,( D 1'^ W p)( D 1^v Z p)- 

To give nonzero contributions, D 2 {ll) must contain three functional derivatives: one functional derivative with respect 
to a f) and two with respect to rj. This gives us three kinds of terms: (i) the argument of fj does not come from the 
same D as the arguments of the two 77, as in D 2 {1 , ) = 61,6262' (ii) the argument of f\ is the same as that of one of the 
two 77, for example D (1 ^ = SiSyS^ (hi) the argument of 77 comes from the same D as the argument of one of the two 
77 but they are not identical, for example D {1 > ) = —616262/. In A f D 2 there are four terms of type (i), eight terms of 
type (ii) and eight terms of type (iii). All the terms of each class give the same contribution. We calculate now the 
contribution of the three kinds of terms 



1 . Terms of kind (i) 

The sum of the terms of type (i) is 




dldl'd2d2't;(l - 1>(2 - 2') ^ ± ±' (6 V 626^6^) (6^,6^6^ Z p ) . 



We compute 

6 V 6 2 62>6xW} = -6x6 v 6262>Wl = ±6r(x,l , ,2,2'), 

and 

616 y 6162' 626yZ p = -6x6262, 6\6y6 y Z p = zG(l, 2, 2', 1, 1', y) 

to obtain 

- % -e j dldl'd2d2y 1 - 1>(2 - 2>(s, 1', 2, 2') ±'G(1, 2, 2', 1, 1', y). 

Terms of kind (ii) 

The sum of the terms of type (ii) is 




dldl / d2d2 / v(l - 1>(2 - 2') ^ ± ±' (6 1 6 2 ^6xW^) (6 V 6 V 61626^ Z p ) . 



We use 

6162*6-2,6^1 = -5z5y5 1 5 2 >W} = ±er(x,2',l,2'), 

and 

6y6y6\6 2 6 2 6yZp = -6i6y626y6 2 6yZ p = zG(l, 1', 2, 1', 2, y). 
Therefore, the fifth term of type (ii) becomes 

-ie f dldl'd2d2'v( 1 - l f )v(2 - 2')r(x, 2', 1, 2') ^ ±'G(1, 1', 2, 1', 2, y). 
^ ±,±' 
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3. Terms of kind ( iii ) 

The sum of the terms of type (iii) is 

/ dldl / d2d2 / ^(l - 1>(2 - 2') J2 ±±f {^2k>5xWl)(8y8 v 8i8 2 >k5yZp). 
J ±,±' 

We have 

Sihh'SxW} = SzSvS&W} = Ter(z,2',l,2), 

and 

5i>5i, 6j6 2 ' S^SyZp = — didyd^y 82' dyZp = iG(l,l' ,2,1' ,2' ,y). 
The fifth term of type (iii) becomes 

-ie I dldl / d2d2X 1 - 1>(2 - 2>(x, 2', 1, 2) ^ ±'G(1, 1', 2, 1', 2', y). 
H. Sixth term of equation (3) 

The sixth term of equation (3) is 

It looks complicated because we have now A'D 3 = (A' D)(A' D)(A' D). The third A' D will have the variables 3, 3' 
and ±" . Again, we need one derivative with respect to fj and two derivatives with respect to 77 in D 3 1>y There are 
24 terms of this kind in A'D 3 , and they all give the same contribution as 816263 6i6i>6if 62' 6263' 63/ 63. Thus, the 
sixth term is 

- / dldl / d2d2 / d3d3 / v(l - 1>(2 - 2>(3 - 3') V ± ±' ±"{6i5 2 6 3 5^W^) ((Mi'M^M^a'Ms^p)- 
2 J ±±/± „ 

We compute 

6162636^ = -5z5i5 2 5 3 W} = ±±' ±"er(z, 1, 2, 3), 

and 

6\6v 6y 6 2 ' 62' 6263' 6^ 6^6 y Z p = 6y 626^6363,616^6^ 6^ 6 y Z p = zG(l', 2, 2' 3, 3', 1, V 2', 3', y). 
The last term of equation (3) is then 

--e I dldl / d2d2 / d3d3X 1 - 1>(2 - 2>(3 - 3')r(x, 1, 2, 3) ^ G(l', 2, 2', 3, 3', 1, 1', 2', 3', j/). 
^ ±±'±" 




+ 





FIG. 1: Diagrammatic representation of equation (10). 



I. The Green function hierarchy 

If we gather all the previous results, we obtain the Green function hierarchy 

G(x,y) = G^{x,y)=G% € ,(x,v)-i fdldl f v(l-l f )J2^G° p (x,l)G(l,l\l\y) 

J ± 

-ice' j dldiy 1 - l')r(x, l 7 , 1, y) ^ ±G(1, 1') + iee' f dldl^(l - l')r(x, 1, 1, 2/) ^ ±G(1 7 , l 7 ) 
J ± J ± 

-ie [ dldl / ^(l-l / )r(x,l,l,l / )^G(l / , 2/ ) 

+ -e / dldl / d2d2 / v(l - 1>(2 - 2>(x, V , 2, 2 7 ) ^ ± 7 G(1, 2, 2 7 , 1, V \ y) 
2 ^ ±,±' 

+e / dldl / d2d2 / v(l - 1>(2 - 2 7 )r(z, 2 7 , 1, 2) ^ ± 7 G(1, l 7 , 2, l 7 , 2 7 , j/). 
^ ±,±' 

-fee 7 / dldl / d2d2 / v(l - 1>(2 - 2>(z, 2, 1, y) ^ G(2 7 , 1, l 7 , 2, 2 7 , l 7 ) 
^ ±±' 

+e / dldl / d2d2 / v(l - 1>(2 - 2')r(x, 2 7 , 1, 2 7 ) ^ ± 7 G(1, l 7 , 2, l 7 , 2, y). 

±,±' 

--e / 4 dldl 7 d2d2 7 d3d3 7 v(l-l 7 )v(2-2 7 )v(3-3 7 )r(x,l,2,3) ^ G(l 7 , 2, 2 7 , 3, 3 7 , 1, l 7 , 2 7 , 3 7 , y). (10) 

^ ±±'±" 

In terms of the Feynman diagrams of nonequilibrium quantum field theory 2 , these terms can be rewritten as in 
figure 1. The order of the terms is the same in equation (10) and in figure 1. 

IV. CONCLUSION 

The hierarchy of Green functions presented here is an important step in the effective calculation of degenerate 
systems. However, the practical implementation of this hierarchy requires an approximation to close it. The GW 
approximation 3 5 is a well-known and powerful way of closing the hierarchy. It has to be adapted to our more 
complicated hierarchy. For the calculation of the optical properties of allochromatic crystals, it will also be necessary 
to adapt the Bethe-Salpeter approach 6 to our setting. 
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Moreover, a functional derivation of the energy with respect to the density matrix provides equations that enable us 
to unify the Green-function formalism and the diagonalization method of many-body theory. This will be presented 
in a forthcoming publication 7 . 
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